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I Introduction 

By a magic square of order n is here meant an arrangement, without 

repeats, of the integers {0, 1, 2, ...n 2 - 1} into the n 2 cells of an n x n square 

^ ■ in a way that the sum of the elements of each row, of each column and of 

each of the two diagonals is the same. Since the sum of 1 + 2 + ... + {n 2 — 1) 

is n 2 (n 2 - 1)/2 and since for a magic square this must be equal to the sum of 

the sum of the n-integers in each of the n-rows (or columns) it follows that 

^ ■ the common sum, cr 2 of these integers must be 

t> 

o: a 2 = — 5 — (l.i) 

O ' Figure I shows two magic squares for n = 4 and n = 5 
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Figure I 



with sum values 30 and 60, respectively. Of particular interest is the magic 
square in lb. Not only do the sums along the rows, the columns and the 
two diagonals have the value 60, but so do the eight broken diagonals that 
are obtained if all partial diagonals are completed by imagining the square 
bent around into a cylinder. Examples of broken diagonals in Figure lb are: 
21, 10, 4, 18, 7; 9, 23, 12, 15, 1; 17, 6, 3, 14, 20; 11, 22, 8, 19, 0. Magic 
nxn squares, such as this one, where the sums along the n-rows and columns, 
the two diagonals and the 2n — 2 broken diagonals all have the same value are 
called pandiagonal or "diabolic" magic squares [9, ch. 10], [11]. To simplify, 
we shall use the term ROW in all of the following to mean any one of a row, 
a column, a diagonal or a broken diagonal. Thus in a pandiagonal magic 
square the sum along any ROW is the same and given by (1.1). 

By a magic pandiagonal cube of order n we shall mean an arrangement 
of the integers 0, 1,2, ...n 3 - 1 on the lattice points of an n x n x n cube 
so that the 3n squares parallel to the faces of the cube as well as the six 
"diagonal" squares which bisect the cube and contain its 4 body-diagonals 
are all pandiagonal magic squares. Thus in each of the 3n + 6 (overlapping) 
squares contained in the cube, the sums along all of the the ROWS — which 
here includes also the files — are the same. It is easily seen that in this case 
the common sum cr 3 must be 

n(n 3 — 1) . 
^ = { 2 } (1-2) 

By extension, we define a four dimensional (4-D) pandiagonal magic 
n x n x n x n cube as one in which the integers 0, 1, 2,...n 4 — 1 are placed, 
without repeats, at the n A lattice sites of a 4-D hypercube, so that the three 
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dimensional (overlapping) n x n x n cubes that can be formed within it, are 
pandiagonal magic cubes. The sum along each ROW of a 4-D cube is easily 
shown to be 

a 4 = n(n 4 - l)/2 (1.3) 



In constructing pandiagonal magic squares, cubes, and hypercubes, we 
follow the idea of Euler [6] — who was concerned exclusively with magic 
squares — and used latin squares (LS) in their construction. As defined 
by Denes and Keedwell [4] and by Laywine and Mullen [9] a latin square 
of order n is an n x n array each of whose lattice points is occupied by 
one of n given symbols, in a way so that no row or column contains any 
one of these symbols more than once. The number of LS's grows rapidly 
with n and for example for n — 10 and 15 this number is ~ 10 36 and 
~ 1.5 x 10 86 , respectively [9, p. 5]. Here we are interested in the relatively 
small subset of the LS's, the pandiagonal latin squares, for which in ad- 
dition each diagonal as well as each of the 2n — 2 broken diagonals — in short 
it's ROWS — also contains each of the n elements precisely once. 

For a general algebraic theory of pandiagonal (diabolic) magic squares 
see the analysis by Rosser and Walker [11]. The books by Andrews [2], 
Kraitchik [5] and Rause Ball and Coxeter [3] contain a more empirical ap- 
proach to magic squares and include various examples and practical rules 
for constructing magic squares. Pasles [10] describes some unusual magic 
squares constructed by Benjamin Franklin. 

Some 5x5 pandiagonal magic squares and some of their important 
features are described by Gardner [7]. In a companion column [8] he details 
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a magic cube of order seven and notes its non-pandiagonal nature. This is 
reenforced by the results of Wynne [12] who studied magic cubes of order 
seven and showed that even if every square of the cube that is parallel to 
a cube-face is pandiagonal, not all of the six, diagonal squares of the cube 
can be pandiagonal. This is consistent with the present analysis, according 
to which only for n > 11 do magic pandiagonal cubes exist. (Similarly we 
find that only for n > 17 can 4-D pandiagonal magic cubes exist.) That 
non-pandiagonal magic cubes of order 7 do exist, however, is established by 
Andrews [2] and by Alspach and Heinrich [1], the latter, incidental to their 
discussion of cubes of order 4m. 

II. Odd-n Pandiagonal Latin Squares and Cubes in 3 and 4 Dimen- 
sions 

A. Pandiagonal Latin Squares 

To set the stage for our discussion of magic pandiagonal cubes in 3 and 
4 dimensions, in this section we collect — and in some instances amplify on— 
properties of pandiagonal latin squares [4], [9]. 

A latin square of odd order n, LS, is an n x n array involving n distinct 
symbols with each symbol appearing once and only once in each row and 
column. We shall invariably take the n symbols to be the integers comprising 
the set S defined by 



S = {0,1, 2, ..,n-l} 



(2.1) 



Also useful will be the set S defined by 



S = {l,2,..,n-1} 



(2.2) 
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A diagonal LS has the additional property that each of its two diagonals 
also contains the n chosen symbols exactly once. Finally a pandiagonal LS 
is a diagonal LS in which, in addition, each of the n symbols appears once 
and only once in each of the 2n — 2 "broken diagonals" , as defined in the 
preceding section. 

Figure II shows 4 LS's of order 3,4,5, and 5, respectively. The 3 x 3 LS 
in part (a) is not a diagonal LS by virtue of the three zeros 
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appearing in the upper left to lower right diagonal. By contrast the 4x4 
LS in part (b) is diagonal but not pandiagonal while the 5 x 5, LS in 
part (c) is pandiagonal. Note that, given a pandiagonal LS, if we per- 
mute its symbols {0, 1,2, ...n — 1} it remains pandiagonal. For example, 



if in lie, we carry out the permutation — > 1; 1 — > 2; 2 ^ 3, 3^4; 
4^0, the resultant LS remains pandiagonal and is given by Figure lid. 
Thus a given n x n pandiagonal LS is the basis for n\ different ones that 
result from the n\ possible permutations of the n-symbols. 
Consider for odd n, the n x n array L with elements Ljj 

Lij = oi\i + a 2 j (mod n) (2.3) 

where i,j run over the elements of S as do the elements themselves and 
a-y and a 2 are non-zero positive integer parameters and thus are elements of 
S in (2.2). When we wish to stress the dependence of L on a\ and a 2 , we 
shall also use the notation L — < cci, a 2 >■ 
We now establish: 

Theorem (2.1): If L — < a\,a 2 > is the array in (2.3) and the greatest 
common divisor with n of each of «i, a 2 , a\ ± a 2 is 1, that is 

(«i,n) = l; i = l,2 (2.4a) 

(ai + a 2 ,n) = 1 (2.46) 

(«i — CK2,ri) = 1 (2-4c) 

then L is a pandiagonal LS for n > 5. (Note that (2.4) cannot be satisfied 
for even n.) 

Proof: Firstly, since for n — 3, the possible a-values are 1 and 2 and since 
these violate (2.4b), it follows that (2.3) is not a pandiagonal LS for n — 3. 
On the other hand, forn = 5, the pairs [cci, a 2 ] = [1,2] and [1,3], for example, 
do satisfy each of (2.4). 
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Secondly to establish first that L in 2.3 for n > 5 is a diagonal LS, consider 
it for fixed j, say, as i ranges over the jth row. As i thus varies from 
to n — 1, the jth row of L varies over the same set by virtue of the hypothesis 
(eti,n) = 1. Similarly for fixed i, as j varies over the ith column, since 
(a 2 ,n) = 1, Lij varies over the same set S. Thus the rows and columns of 
L satisfy the condition that L be a latin square. Along the diagonal from 
(i, j) = (0, 0) to (n — 1, n — 1), i = j, so that along here = (a± + 012)1 
(mod n). Thus again as i varies over the set S the n diagonal elements of L 
must be some permutation of S since (ct\ + a2,n) = 1 according to (2.4b). 
Finally along the other the diagonal from (i, j) = (0,n — 1) to (n — 1,0), 
i + j — n — 1, so that = (a\ — 012)1 + «2(n — 1) mod n which has the same 
property by virtue of (2.4c), (a^ — a 2 , n) = 1. Thus we have established that 
L in (2.3) under the constraints in (2.4) is a diagonal LS; it remains only to 
establish that it is also pandiagonal. 

To this end, consider, for example, the "split" diagonal just above — and 
parallel to — the lower left to upper right diagonal of L and its appendage 
in the lower right hand corner, (i, j) = (n — 1,0). (We assume i increases 
from to the right and j increases upward from 0.) Its entries are defined 
by j = i + 1, i = 0, 1, 2, ...n — 2 and (i, j) — (n — 1, 0). Substituting j — i + 1 
into (2.3) we find 

L i>i+ i = («! + a 2 )i + a 2 (mod n) 

If we now let i run over the complete set S, it is easy to see that since 
Li n = Li (mod n) we obtain both parts of the broken diagonal and that no 
two of these n-elements are the same since they are simply a permutation of 
the elements of S by virtue of (2.4b). This argument is easily repeated for 
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all broken diagonals and we conclude that each broken diagonal consists of 
a permutation of the elements of S. We have thus established the theorem. 
Remark 2.1: If a±, a 2 satisfy each of (2.4), then so do ka±, ka 2 for any keS 
for which (k, n) = 1. 

Remark 2.2: Reference to (2.3) shows that if L =< ai,a 2 > is a pandi- 
agonal LS, then so is < ka±,ka 2 > for any positive integer keS provided 
(k, n) = 1. For according to (2.3), < ka±,ka 2 >, is simply < ai,a 2 > with 
its elements permuted in some way. 

Remark 2.3: If we add an integer xeS to the right hand side of (2.3), we 
obtain < a±, a 2 > with its elements permuted in some way. 

An important notion relating to LS's is that of the orthogonality [9, ch. 
2]. Two n x n, LS's are said to be orthogonal, if when they are superposed, 
none of the n 2 ordered pairs of elements that result occurs more than once. 
Thus if L^f and L^f eS are any corresponding elements of the two LS's, 
and then they are orthogonal if and only if the ordered pairs [L^\ L^f] 
and [I/^g, L^p] differ (mod n) for any fixed values of i, jeS but for all choices 
of a, (3, eS. 

Remark 2.4 If two LS's are orthogonal, they remain so when the elements 
of either or both undergo arbitrary permutations. 
We now establish: 

Theorem (2.2) Let =< ai,a 2 > and L (2) = < /3i,/3 2 > be two pan- 
diagonal LS's with both pairs a±,a 2 and (3 2 each satisfying (2.4). If the 
determinant d 2 defined by 



do 



ol\ a 2 
Pi fa 



(2.5) 



8 



is relatively prime to n, i.e. 



(d 2 ,n) = l (2.6) 
then L^> and are orthogonal pandiagonal LS's. 

Proof: Suppose for i,jeS there is a second pair k, £eS for which the ordered 
pairs [LJ^, Lij^] = [L^,L^] mod (n). That is, suppose 

ct\i + a^j = + a 2 £ 

(mod n) 

fci + fcj = Pik + foi 
which for convenience we express in matrix notation 



ol\ a 2 ^ 



A A 



2 



(J " «) 



= (mod n) 



Now since by hypothesis (d2,n) = 1, so that in particular d 2 ^ (mod n), 
we can multiply this relation by 

d ( "i "2 _/ A -«2 \ 
\ A A / V -A ai / 



d 2 =0 mod n 



so that it becomes 



V - Q 

Finally, since (d 2 ,n) = 1 it follows that % — k and j = i. Thus and 
are orthogonal pandiagonal LS's. 

The orthogonality criterion of the two LS's in (2.6) is very convenient and 
as we shall see is extendable to higher dimensions. 

Let us consider the question as to the number of distinct pairs [0:1,0:2] 
there are for given n for which L =< a±,a 2 > is a pandiagonal LS. To 
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simplify let us assume in the following that n is a prime p. For a given prime 
p > 5, consider the p — 3, a 2 ], pairs 

Kc^] = [1,2], [1,3], [1,4],. ..[l,p- 2] (2.7) 

Obviously each of these satisfies (2.4) and no other pair for which a± — 1 
does so. 

Theorem (2.3): The pandiagonal latin squares associated with the [0:1,0:2] 
pairs in (2.7) are: 

(1) pairwise mutually orthogonal 

(2) any other pandiagonal LS can be obtained from one associated with 

< 1, £ >, £ = 2, 3, ...p — 2, in (2.7) by a permutation of symbols. 

Proof: 

(1) Consider the two pandiagonal LS's < 1, £ >, < 1, m >; I, m — 2, 3, ...p — 

2, (£ 7^ m). According to (2.5) the associated determinant d 2 is 

1 £ 

g?2 = = m — £ 

1 m. 

and obviously satisfies (2.6) since we assume m 7^ £. 

(2) Let x, y be any two unequal positive integers eS = {1,2, ...p — 1} that 

satisfy each of (2.4) so that < x, y > is a pandiagonal LS. If x' 1 is the 
inverse of x (mod p) then making use of (2.3) and Remark (2.2), we 
find that < 1, x~ l y > is also pandiagonal and is obtained from < x, y > 
by a permutation of its symbols. Finally since x~ x y (mod p) must be 
one of (2, 3, ...p — 2) in (2.7) it follows from Remarks 2.1 and 2.2 that 
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< x,y > can be obtained from < l,x 1 y > by a permutation of its 
symbols. 

For example, for p = 7, successive multiplication (mod p) by use of k = 
2, 3, ...6, leads to < 1, 2 > -> < 2, 4 >, < 3, 6 >, < 4, 1 >, < 5, 3 >, < 6, 5 >. 
Similarly < 1, 3 > — > < 2, 6 >, < 3, 2 >, < 4, 5 >, < 5, 1 >, < 6, 4 >; and 
< 1,4 > ->< 2,1 >, < 3,5 >, < 4,2 >, < 5,6 >, < 6,3 >. 
Remark 2.5: If we allow for permutation of symbols then all pandiagonal 
LS's in (2.3) can be obtained, for given p by use only of the [a±, a 2 ] pairs in 
(2.7). 

Finally, as shown for pandiagonal magic squares, by Ball and Coxeter [3, 
p. 203] and by Martin Gardner [7] for a 5 x 5 pandiagonal magic square, we 
find that pandiagonal LS's have an analogous unusual property. If in an n x n 
pandiagonal LS we move the left hand column so that it becomes the right 
hand column (or vice versa) or similarly move the top row to the bottom, 
the resultant array is again a pandiagonal LS. This has the consequence that 
if we "tile" the plane with a given n x n pandiagonal LS we can outline any 
n x n square on this infinite pattern and obtain a pandiagonal LS. 

The underlying result is contained in: 
Theorem 2.4: If the left column of a pandiagonal latin square is moved so 
it becomes the right column, the resulting LS is a pandiagonal LS obtained 
from the original one by a permutation of symbols. 

Proof: Let L =< a±,a 2 > be the original pandiagonal LS and define a 
second one L' with elements given by 

L[ ■ = Lij + ai = («i + l)i + a 2 j (mod n) 
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It follows from (2.3) that V is simply L with its elements permuted. Further 
we have 

L 'oj = L if L 'ij = L 2j] - L n-i,j = L of {mod n); j = 0, 1, ...n - 1 

which shows that L' is simply L with its left column moved so it becomes 
the right column. 

For example, if the left hand column of the LS in Figure lie, which inci- 
dentally is simply < 1,3 >, (with i increasing to the right and j increasing 
upwards) is moved to the right side, the original LS, but with permuted sym- 
bols, in lid results! Of considerable interest perhaps is that a form of this 
property, as will be shown below, has an analogue in higher dimensions. 
B. Pandiagonal Latin Cubes 

We define a pandiagonal latin cube asannxnxn cube each of whose 
n 3 lattice points contains one of the members of S = {0, 1,2, ...n — 1} and 
in a way so that each of its 3n + 6 constituent squares is a pandiagonal latin 
square. Recall in this connection that an n x n x n cube has 3n squares 
parallel to a cube face plus 6 "diagonal" squares which contain its 4 body 
diagonals. Note that, as here defined, in a pandiagonal latin cube each row, 
column, file, diagonal and broken diagonal of each of its squares, i.e., its 
ROWS, contain each element of S once and only once. 

In an obvious generalization of (2.3) to three dimensions, we define an 
n x n x n array C by the formula 

djk = ol\% + «2j + 0:3k (mod n) (2.8) 

where each element of e5,{0,l,...n-l}. Here i,j,k are integer variables 
each running over {0,1,. ..n — 1}, and ai,a 2 ,a 3 are elements of 

12 



S = {1,2, ...n — 1}. We shall use the notation C = < 01,02,03 > when 
we wish to focus on the dependence of C on 01, 02, 03. We now establish the 
following: 

Theorem (2.5): The cube defined in (2.8) is a pandiagonal latin cube — 
in that each of its constituent 3n + 6 squares is pandiagonal — provided 



01,02,03 satisfy the constraints: 

(a e ,n) = l; £ = 1,2,3 (2.9a) 

(a t + a e >,n) = l; £, V = 1, 2, 3; t ^ £' (2.96) 

(a t -a e ,n) = l; £, t = 1, 2, 3; t ^ £' (2.9c) 

(A,n) = l; A = a 1 + a 2 + a 3 (2.9d) 

(A-2a e ,n) = l; £=1,2,3 (2.9e) 



Proof: Consider first the 3n squares parallel to the faces of the cube. For 
fixed k, say, < k < n — 1, consider the square of that is parallel 
to the i — j plane as i and j vary over the elements of S. Then has 
essentially the same structure as does in (2.3) since the constant a 3 k is of 
no consequence. Making use of the restrictions on the o's in (2.9 a,b,c) and 
comparing with those in (2.4 a, b, c) we conclude that for any fixed k, the 
square is a pandiagonal latin square. Repeating this argument for fixed 
i, with j and k variable and for fixed j with i and k variable we conclude 
that all 3n squares in the cube defined by (2.8) that are parallel to a cube 
face are pandiagonal latin squares. 
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With regard to the six squares not parallel to a cube face we proceed 
as follows. For that diagonal square, with vertices at (i,j,k) = (0,0,0), 
(0, 0, n — 1), (n — 1, n — 1, n — 1), (n — 1, n — 1, 0), we have i = j so that (2.8) 
becomes 

Cuk = («i + a 2 )i + a 3 k (mod n) . 

But this is again of the form in (2.3) if we make the replacements in (2.4a,b,c) 
«i — > cti + a 2 ; «2 — ► «3- The first of (2.4) is satisfied because of (2.9a,b) 
while (2.4b) and (2.4c) become respectively (a\ + a 2 + «3,n) = 1 and (a± + 
a 2 — a 3 ,n) = 1 and these are the same constraints as in (2.9d) and (2.9e), 
respectively. Thus Cuk, the given square, is a pandiagonal latin square. 

Similarly for the diagonal square perpendicular to Cuk whose vertices 
have the coordinates k) = (n — 1, 0, 0), (0, n — 1, 0), (0, n — 1, n — 1), 
(n — 1, 0, n — 1), we have % + j = n — 1 and 

Ci^ n -i-i),k = ("i - OL 2 )i + a 3 k + a 2 (n - 1) (mod n) 

which on comparison with (2.3) and (2,4) with the replacements a 1 — > a\ — 
a 2 ; a 2 — ■> a 3 leads to the conditions (cti — a 2 + 0:3,71) = 1, and (o;i — a 2 — 
0:3, n) = 1 The first of these is the same as (2.9e) with 1 = 1 and the second 
is the same as (2.9e) with I = 1 if we make use of the fact that (-x, y) = 1 
is equivalent to (x,y) = 1. 

A similar argument shows that for the constraints in (2.9) the remaining 
four "diagonal squares" of the cube are also pandiagonal latin squares. The 
theorem is thus established. 

Just as for the 2-D case, we require a 3-D analogue of orthogonality of latin 
cubes. For our purposes, we shall say that three latin cubes are orthogonal 
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if when they are superposed none of the n 3 ordered triplets of elements that 
result occurs more than once. There are other definitions of orthogonal latin 
cubes, [9, ch. 3], but for purposes of producing magic cubes this definition 
is essential. 

We now establish: 
Theorem (2.6): Consider the 3 pandiagonal latin cubes 



(mod n); 



9=1,2,3 



(2.10) 



where a pq (p,q = 1,2,3) are elements of S and for each value for q satisfy 
the conditions in (2.9) and let d% be the determinant d% = \a pq \. Then if d^ 
is relatively prime to n, that is 



(d 3 ,n) = 1 

then the three cubes C^ q \ q = 1, 2, 3 are orthogonal. 



(2.11) 



Proof: Suppose on the contrary that for a given (i,j,k) there existed an 
integer triplet (u,v,w) each eS not equal to (i,j,k) for which C^ k = C q vw 
(mod n) for each q. Then, as for the 2-D case, we could express this in matrix 
notation 

/ 

an 0£ 2 i a 31 

au «22 «32 3 — V 

k — w 

Since by (2.11), the determinant d 3 is relatively prime to n, so that in par- 
ticular d 3 ^ mod n, we may multiply both sides of (2.12) by d 3 times the 



\ ( i-u \ 







mod n 



(2.12) 
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inverse of the matrix on the left. The result is 




V 



k — w 



Finally since (d 3 ,n) = 1 it follows that i = u; j = v; k = w and the theorem 
is proved. 

An empirical study of the constraints in (2.9) shows that no triplets 
[ai, a 2 , ct 3 ] satisfying (2.9) exists for n = 5,7,9. (The latter clearly since 
n would be divisible by 3.) Such an analysis is most easily carried out 
by recognizing that — as in the 2-D case — without loss of generality we can 
take ot\ = 1 and consider simply the triplets [1, 0:3]. For n — 7, it is 
easily confirmed that no values, such [1, 2, 3], [1, 2, 4], [1, 2, 5] would satisfy 
all of (2.9). For n — 11 however we find, among others, the possibilities 
[1,2,4], [1,2,5], [1,2,6], [1,2,7], [1,5,8], [1,6,8] as well as these with a 2 and 
0:3 interchanged. As an example of orthogonal pandiagonal cubes we note that 
for appropriate integers, £,m,p for the three LS's < 1, 2, £ >, < 1, m, 2 > 
and < 1, 2,p >, for n = 11, d 3 is given by 



and will for values of £, m, p with £ 7^ p and 4 < £, m, p < 9 lead to (0^3, n) — 1. 
Figure III shows two planar sections through the cube < 1, 2, 7 > for n — 11 : 
(a) corresponds to the square k = 2 in (2.8) and (b) to the diagonal square 



1 2 £ 



1 m 2 



(m-2)(p-£), 



1 2 p 



16 



i = j. 

Note, however, that for the triplet < 1,2,£ >, < l,2,m >, < l,2,p >, 
d 3 = so that these three do not constitute orthogonal cubes as we have 
defined them even though each cube is itself pandiagonal. 

As for the analogous 2-D case (Theorem 2.4) we can easily establish the 
fact that if we move, say, a face of a pandiagonal latin cube to its opposite 
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Figure III 



side, the resultant cube remains a pandiagonal latin cube. To see this, con- 
sider C = < ai,a 2 ,a 3 > in (2.8) and define C with elements given by 
C' ijk = C^k + on so that 
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C' ijk = ai(i + 1) + a 2 j + a 3 k (mod n). 

Obviously, 

c 'ojk = c ijk,C l ljk = C 2jk -,...C l n ^ ljk = C 0tjk (modn); < j, k < n - 1 

Thus C" is obtained from C by transporting its i — face to i — n — 1. This 
implies that if we "tile" all of 3-D space with a given pandiagonal latin cube 
any n x n x n cube selected out of this infinite array will be a pandiagonal 
latin cube but with its elements permuted in some way. 
C. Four dimensional pandiagonal latin cubes 

A 4-D pandiagonal latin cube is an arrangement of the integers 
{0, 1,2, ...n — 1} among the cells of an n x n x n x n cube in a way so 
that each of its 4n + 12 constituent 3-D cubes is a pandiagonal latin cube. 

In an obvious generalization of (2.8), we define an n x n x n x n array of 
integers by the formula 

Hijkf = a\i + a 2 j + a s k + a^£ (modn) (2-13) 

where each of H^m is e S = {0, 1, 2, n — 1} and where k, £ are integer 
variables each with the range 0, 1, 2, n—\. The four quantities a±, a 2 , 0:3, 0:4 
are integer parameters e S — {1,2, ...n — 1}. As before we use the notation 
< cti, a 2 , a 3 , «4 > when the dependence of H on the a's is of interest. We 
now establish the following: 

Theorem (2.7): The hypercube defined in (2.13) is a 4-D pandiagonal 
latin cube provided the four integer parameters cti, a 2 , a 3 , a 4 satisfy the con- 
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straints 



(a m ,n) = 1 



m=l,2,3,4 (2.14a) 

771,771' = 1,2, 3,4; m^m' (2.146) 

m,m' = 1,2, 3,4; m ^ w! (2.14c) 

B = ai + a 2 + « 3 + a 4 ;m = 1,2,3,4 (2.14d) 




(2.14e) 



( J B-2a m ,n) = 1 

(B - aw - 2a m ,n) = 1 

(5 - 2aw - 2a m , n) = 1 



m = l,2,3,4 (2.14/) 
m,m! = 1,2, 3,4; m ^ m' (2.14#) 
m,m' = 1,2, 3,4; m ^ m' (2.14/i) 



Proof: Consider first the 3-D cube that results for fixed £ from (2.13) as i, j, k 
vary over S. It will be a pandiagonal latin cube provided the constraints 
in (2.9a)-(2.9e) are satisfied. Now (2.9a)-(2.9e) are contained within (by 
appropriate choice of the subscripts) (2.14a), (2.14b), (2.14c), (2.14e) and 
(2.14g), respectively. Thus since the added constant plays no role the 
cube is a pandiagonal latin cube. Similarly, the cubes that result for fixed i, 
as j, k, £ vary and for fixed j as i, k, £ vary and for fixed k as i, j, £ vary over 
S are all pandiagonal latin cubes. There are altogether 4n pandiagonal 3-D 
latin cubes of this type. 

Similarly the "diagonal" 3-D cube that results for i = j, as i,k,£ vary 
over S has the form 



On comparison with (2.8) and (2.9a) - (2.9e) we see that this is also a pandi- 
agonal latin cube since we assumed (2.14a), (2.14b) (2.14c), (2.14d), (2.14e), 
(2.14g) to be satisfied . And similarly for the other five pairs: % = k,i = £, 



H iiki = («i + a 2 )i + 0.3k + a^£ {mod n) 



j = k,j=£,k = £. 
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Finally for the "diagonal" latin 3-D cube that results for % + j = n — 1, 

as i, k and £ vary over S 

Ei, n -\-i,k,t = («i - a 2 )i + a 3 k + a 4 £ - a 2 (n - 1) (mod n) 

which on comparison with (2.9a) - (2.9c) is a pandiagonal latin cube by virtue 
of (2.14a), (2.14b), (2.14c), (2.14f), (2.14g), (2.14h). And similarly for the 
remaining five cubes i + k — n — l,i + £ — n — 1, j + k — n — 1, j + £ — n — 1 
and k + £ — n — 1. This proves the theorem. 

Turning to the question of orthogonality, we define four, 4-D pandiagonal 
cubes to be orthogonal, if when they are superposed no two of the n 4 ordered 
quartets of elements that result are the same. As in the lower dimensional 
cases there are other definitions of this orthogonality [9] , but for our purposes 
this one is essential. 

By analogy to Theorems (2.2) and (2.6) we have: 
Theorem (2.8): Consider the four, 4-D pandiagonal latin cubes 

H-f ke = Oii q i + a 2q j + d 3q k + a iq £ (mod n); q — 1, 2, 3, 4 (2.15) 

where a pq (p, q = 1,2, 3, 4) are elements of S which for each value of q satisfy 
the conditions in (2.14) and the variables i,j,k,£ range over 
S = {0, 1,2, ...n — 1}. Further, let d 4 = \a pq \ be the determinant of the 
a's. Then if d 4 is relatively prime to n 

(d A ,n) = l (2.16) 

then the four 4-D cubes in (2.15) are orthogonal. 

Proof: Suppose to the contrary there existed a quartet of integers (u, v, w, x) 
each e S, not equal to any given quartet k,£) for which H^l e = H^ wx 
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(mod n) for each q. Then just as in deriving (2.12) we would find 



an Ol 2 l «31 «41 

CK12 «22 «32 «42 

"13 "23 "33 "43 

y «14 «24 «34 «44 / 



\ 







mod n 



(2.17) 



J - V 
k — w 
\i-x ) 

and conclude following essentially the same steps as before that k, £ must 
be equal to u, v, w, x respectively. Thus concluding the proof of orthogonality. 

It is easily confirmed by enumerating the various possibilities that only 
for n > 17 is it possible to find integers [aii, ol^ a^, 0:4] among {1, 2, ...n — 1} 
for which (2.14) can be satisfied. In particular for n — 17, possible 4-D 
latin hypercubes are given by < 1,2,4,8 >, < 1,2,4,9 >, < 1,2,13,8 >, 

< 1, 2, 13, 9 > as is readily confirmed. Thus a possible form for the determi- 
nant of the matrix in (2.17) is 

12 4 8 
12 4 9 
12 8 4 
14 9 2 

and since (-8, 17) = 1, the four 4-D pandiagonal latin cubes < 1,2,4,8 >, 

< 1, 2, 4, 9 >, < 1, 2, 8, 4 >, < 1, 4, 9, 2 > constitute an orthogonal set of such 

hypercubes, for n — 17. 

Figure IV, shows a planar section through the hypercube < 1,2,4,9 > 
corresponding to % = 2, j + k = 16 in (2.13). 
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Figure IV 



Just as for the 2 and 3 dimensional cases, if H — < ati, «3, «4 > is a 
4-D pandiagonal latin cube, then so is H' with elements given by 

H'ijkt = H ijk e + ol\ = («i + l)i + a 2 j + a 3 k + a 4 £ (mod n) . 

Again we can imagine "tiling" all of 4-D space with H and be assured that 
any nxnxnxn subcube in this space will be a 4-D pandiagonal latin cube, 
with its elements a permutation of the original elements of H . 
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III. Magic, Pandiagonal Squares and Cubes in Two and Three Dimen- 
sions 

With the results of the preceding section available, it is now straightfor- 
ward [3], [6] to generate magic pandiagonal cubes in three and four dimen- 
sions. To set the stage we first illustrate the matter in two dimensions. 
A. Magic Pandiagonal Squares 

Recall that a magic pandiagonal square of order n is an arrangement, 
without repeats, of the integers (0, l,2,...n 2 — 1) on the n 2 lattice points 
of an n x n array so that the sum of the elements in each of the n— rows, 
n— columns, n— diagonals (including the n — 2 broken diagonals) has the same 
value. (As above let us use the generic ROW to represent any one of these 
rows, columns, diagonals, etc. ) This common sum of the ROWS has been 
given in (1.1). 

Theorem 3.1 [9, p. 178] Let = < a u a 2 >, L (2) =< Pi, (3 2 > be two 
orthogonal pandiagonal latin squares that separately satisfy the conditions 
of theorems (2.1) and (2.2) and define an n x n array with elements 

Mlf by 

M< 2 > = nL« + (3.1) 
Then is a magic pandiagonal square. 

Proof: Since the elements of and range over S — {0, 1, 2, ...n — 1} it 
follows from (3.1) that each of the elements must be one of the integers 
0, 1,2, ...n 2 - 1. Further, since the sum of the elements in each ROW of the 
pandiagonal and L^ 2 \ isO + l + 2 + ... + n — 1 = n(n — l)/2 it follows 
that the sum of each of the 4n ROWS of M (2) is 

n[n(n - l)/2] + n(n - l)/2 = n(n 2 - l)/2 = a 2 



24 



with a 2 defined in (1.1). Finally since and are orthogonal, no two 
elements of can be the same and the theorem is established. 

Since and L^> are pandiagonal LS's for which we know n > 5, only 
for n > 5 can in (3.1) be a magic pandiagonal square. Further if we 
tile the plane with any n x n subsquare will also be a pandiagonal 

magic square since L^> and have this same tiling property [7]. (See the 
discussion at the end of II A.) 

Since the elements of the orthogonal pandiagonal latin squares and 
may be permuted among themselves without changing their essential 
properties, it follows that from a given pandiagonal LS we may generate n\ 
different versions. This leads to the number N 2 of pandiagonal magic squares 
obtainable by our method to be 

N 2 = £{n)(n\f (3.2) 

where £(n) is a low order polynomial in n. For n — 5, since the only inde- 
pendent LS's are < 1,2 > and < 1,3 > (Theorem 2.3), and since we may 
interchange their roles in (3.1) £(5) = 2 and we obtain consistent with the 
results of Rosser and Walker [11] 

N 2 = 2880 

Clearly because of the (n!) 2 factor N 2 rises very rapidly with n. For n — 11 
for example, N 2 is ~ 9.0 x 10 16 . 
B. Pandiagonal Magic Cubes 

A magic pandiagonal nxnx n cube is an arrangement, without repeats, 
of the integers 0, 1, 2, ...n 3 — 1 onto the n 3 lattice sites of a cube so that the 
sums along each ROW (i.e. along each of the n 2 -rows, n 2 -columns, n 2 -files 
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and the n(3n + 6) diagonals (including the broken diagonals) in each of its 
(3n + 6) squares are the same. The common sum along the ROWS is given 
by 0-3 in (1.2). 

Theorem (3.2): Let = < a± q , a 2g , a 3q >, q — 1,2,3, be three orthogo- 
nal pandiagonal nxnx n latin cubes that satisfy the conditions of theorems 
(2. 5), (2. 6) and define an n x n x n array with elements by 

Then Af < 3 > is a magic pandiagonal cube. 

Proof: Since the elements of C^,C^ and range over 

S = {0, 1, 2, ...n — 1} it follows from (3.3) that each element of must be 
one of the integers {0, 1, 2, ...n 3 — 1}. Consider now any ROW of M^. The 
sum of the elements of this ROW is, according to (3.3), given by 

n 2 [n{n- l)/2] + n[n(n - l)/2] +n(n- l)/2 = n(n 3 - l)/2 

and yields the value <r 3 in (1.2). Thus the sum of the elements in each ROW 
of the cube is the same and since and are orthogonal, it follows 

from (3.3) that no two elements of can be the same. The theorem is 
thus established. 

According to the discussion in IIB, since the C's are pandiagonal latin 
cubes, for which we know n > 11, will exist only for n > 11. 

As for the two dimensional case, if we allow for the interchange in 
of aig, a 2q , a 3q (q = 1,2,3) with each other, and of permuting the symbols 
in each of and independently we can conclude that N 3 , the 

number of cubes in (3.3) is given by 

N 3 = £ 3 (n)(n\) 3 (3.4) 
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where £3 is an appropriate polynomial in n. For n = 11, 13 and 17, (n!) 3 
assumes the approximate values 6.4 x 10 22 , 2.4 x 10 29 and 4.5 x 10 43 respec- 
tively. Because of this rapid rise of N 3 with n, we can anticipate that the 
factor £3(71) will not affect this variation, qualitatively. 

As for the analogous 2-D case, it follows that by virtue of the "tiling" 
properties of pandiagonal latin cubes (see IIB), we can also tile 3-D space 
with any magic pandiagonal cube, and be assured that any n x n x n cube 
selected out of this infinite array will also be a magic pandiagonal cube. It 
will differ from the original cube in that the elements of its underlying latin 
cubes will have been permuted. 
Four Dimensional Pandiagonal Magic Cubes 

By analogy to the above, we define a magic four dimensional pandiagonal 
cube, as an arrangement, without repeats, of the integers 0, 1,2, ...n 4 — 1 
among the n 4 lattice sites of an n x n x n x n cube so that the sum of the 
elements in each ROW of the 4-D cube, has the same value 04 as given in 
(1.3). 

Theorem (3.3): Let = < a± q , a 2q , «3 g , a iq >, q — 1,2,3,4 be four 
orthogonal pandiagonal n x n x n x n latin cubes that satisfy the conditions 
of theorems (2.6) and (2.7) and define an n x n x n x n array with 
elements M.\f u by 

M< 4 ' = n 3 Q« + n 2 Q^ + n Q& + qW (3.5) 

Then is a magic, pandiagonal four dimensional cube. 
Proof: Since the elements of each of Q (l \Q (2 \Q^ and range over 
S = {0, 1, 2, ...n — 1}, it follows from (3.5) that the elements of M^ M assume 
values from the set {0, 1, 2, ...n 4 — 1}. Further, since the Q's are pandiagonal 
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it follows that the sum of the elements in any ROW of is given according 
to (3.5) by 

n 3 [n(n - l)/2] + n 2 [n(n - l)/2] + n[n(n - l)/2] + n(n - l)/2 = n(n 4 - l)/2 

and this is <r 4 in (1.3). Thus the sum of the elements of any ROW of 
is given by (1.3). Finally because of the assumed orthogonality of the Q's, 
it follows that no two elements of are the same and is a magic 
pandiagonal 4-D cube without repeats. The theorem is established. 

According to the discussion in IIC, since the Q's are 4-D pandiagonal 
latin cubes, for which we found n > 17, it follows that will exist only 
for n > 17. 

The number iV 4 of different M^'s in (3.5) can be estimated as above to 
vary for large values of n as 

N 4 (n!) 4 

so that for n = 17, 19 and 23, N4 assumes the approximate and rapidly 
growing values of n of 1.6 x 10 58 , 2.19 x 10 68 and 4.5 x 10 89 respectively. 

It is also possible, using the tiling properties of the Q's to "tile" 
throughout four dimensional space and obtain a pandiagonal magic 4-D cube 
by selecting any n x n x n x n cube in this space. Such a cube will be the 
same as the original cube but with the elements of its underlying latin cubes 
permuted. 

Evidentally, these arguments are extendable to dimensions higher than 4, 
but the resulting constraints on the a parameters, analogous to those in (2.9) 
in three dimensions and those in (2.14) in four, can be expected to become 
increasingly complex. 
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